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the sine and the cosine of twice the CKM Unitarity Triangle angle /? using a time- 



We present a feasibihty study of a previously outlined method to measure both 



dependent Dalitz plot analysis of B Dtt^tt decays, where the neutral D meson 



in 



is reconstructed in decays to CP eigenstates. We show that this method can be 



used at the B factories to make a measurement of cos(2/3) that is competitive with. 



or more precise than, other techniques using different quark-level transitions, while 



sin(2/3) can be measured to a better precision than any existing measurement using 




b — > cud transitions. Furthermore, this technique has great potential to be employed 
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at LHCb. 
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I. INTRODUCTION 



Precise determinations of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) 

m 

matrix are important to check the consistency of the Standard Model and to search for new 
physics. Measurements of the — mixing phase determined from B^ J/tpK^ (and 
similar) decays at the B factories give sin(2/3) = 0.680 ±0.025 , where 3 is one of the 

angles of the CKM Unitarity Triangle (for an introduction, see, for example [a, [tI, l8|). These 
results confirm the Kobayashi-Maskawa mechanism as the origin of CP violation within 
the Standard Model. Nonetheless, the effects of "new physics" are expected to be seen as 
non- negligible corrections to the Standard Model, particularly if new particles are present at 



energies as low as the TeV scale. 



in several recent reviews 



le purp ose of flavour physics in the LHC 3rcl, clS discussed 
[isl , is to constrain the new physics parameter 



/scale. 1 ne purp ose o 
flQ,[ll|,H[l3[,Q,[l 



space and - once observed - to measure the couplings of the new physics particles. 

An important part of this programme is the precise measurement of the angles of the 
CKM Unitarity Triangle using different quark-level transitions. For example, comparisons 
of the value of sin(2/3) measured in b ^ ccs transitions (such as B^ J/ipK^) with those 
obtained in penguin (loop) dominated h sss transitions (such as B^ — »• (j)K^) or b cud 
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19(1 . In order 



transitions {e.g. B^ Dtx^) can probe for new physics effects 
to achieve the best possible precision, as well as to remove ambiguities in the results, it is 
important to use channels that can measure cos(2/3) as well as sin(2/3). This can be achieved, 
in general, by using any final state that contains interfering amplitudes. 

In this paper, we present the results of a study of the feasibility of measuring cos(2/3) 
from B^ DiT^TT^ decays, where the neutral D meson is reconstructed in decays to CP 
eigenstates. This method has previously been described in outline [2^. In this work, we 
significantly extend the earlier study, taking advantage of results from the B factories that 
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We estimate the sensi- 



provide information on the composition of the Dalitz plot 
tivity that can be achieved with the BaBar dataset, and comment on the potential of the 
LHCb experiment. 

The paper is organised as follows. In Section [III we briefly review alternative approaches 
to measure cos(2/3). The main body of the paper is in Sections IIIII and IIV| in which we 
give a description of the method and discuss the results of our feasibility study. Finally, we 
present our conclusions. 
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II. REVIEW 



The weak phase 2/3 can be probed through mixing-induced CP violation effects in B 
decays mediated by a number of different quark-level transitions. This enables a powerful test 
of the Standard Model, since models of physics beyond the Standard Model that introduce 
new particles at the TeV scale can produce effects that differ between different decay modes. 
In this section we briefly review the various techniques that have been suggested, in order to 
illustrate the need for methods that can be used to provide precise measurements at LHCb 
and future experiments. 

b CCS 

Since the most experimentally precise measurement of sin(2/3) is made using 5° J / ipK^ 
decays, one might expect that a similar channel can be used to extract also cos(2/3). Indeed, 
methods based on — » J / ipK^ decays using the subsequent evolution of the neutral kaon 



have been proposed 
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25| . but appear experimentally challenging. 
The first experimental rneasurement of cos(2/5) used the decay mode 5° J/tpK*{892), 

This method relies on the interference between CP-even and 



with K*{892) i^Ovr" M 



CP-odd helicity states [27|, l28|, l29|. A residual ambiguity due to the unknown sign of 
the strong phase difference can be resolved u sing input from Kn scattering 30|] or from 



theory 



3l|. The most recent measurements 26|, |32| prefer cos(2/5) > but with large 



uncertainty. It will be difficult for LHCb to improve on these measurements since it is 
necessary to measure accurately the momentum of the neutral pion in the final state. 

Other methods using doubly-charmed final states D^*^^ D^*^~ Kg have been proposed. In 
principle, time-dependent amplitude analyses of these states would yield information on 



the weak phase 
po ^ D+D~Ks 



see the discussion in the next section). This is simplest for the decay 
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where all final state particles are pseudoscalars so that the Dalitz 
plot gives a complete description of the phase space. The P factory statistics have not yet 
enabled this analysis, and due to the high multiplicity of charged tracks in the final state it 
may be difficult to study at LHCb. If either or both charmed mesons are reconstructed as 
D* there are additional degrees of freedom that further complicate the amplitude analysis. 
Fortunately, by integrating over regions of the P° D*^D*~ Kg phase space, some sim- 



plifications are possible, but at the price of considerable theoretical 



some input from theory, the current experimental measurements 35|, l36| prefer cos(2/5) > 



uncertainty 



3|. With 
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b — » cud 

The possibility to measure cos(2/5) from Dcpir^'n'' decays [20| is the subject of 

this paper. A similar analysis using D^pn^ 7r~ , with D'^ D±7r^ or -D=p7 [s^] 

is possible in principle, but requires a more involved amplitude analysis including the D* 
decay angles. 

Another method to measure cos(2/3) in 6 — cud transitions, using B^ Dh^ (with hP 
being a light neutral meson such as a tt") with time-dependent Dalitz plot analysis of the 



subsequent neutral D meson decay to KgTi'^n has been proposed [38| and implemented [39|, 



401 ]. To analyze the latter decay channel, the B decay vertex position must be determined 



from the pion tracks that originate from the D decay. The lack of primary particles from 
the B vertex, together with the necessity to reconstruct the neutral meson, will make this 
analysis difficult to carry out in the hadronic environment of LHCb. 

Another interesting approach could be to carry out a simultaneous analysis of the B 
and D meson decay Dalitz plots in the B^ Dn^n^ , D Kgir^ir' decay chain, thereby 
combining the method of Ref. 20] and this paper with that of Ref. js^]- Large statistics 
could be available for such an analysis, since the relevant branching fractions for both B 
and D decays are reasonably high. Although the complete four-dimensional amplitude 
analysis would be quite complicated, it may be possible to select regions of the phase space 
where simplifications are possible (for example, selecting the Dp^ dominated region of the 
B decay phase space). In order to reach high precision, however, a complete modelling of 
the amplitude is likely to be necessary. 

b ccd 

Several possibilities to measure cos(2/5) from b — ^ ccd transitions have been discussed in 
the literature, but all are experimentally challenging and none have yet been implemented. 
For example, the interference between B^ D**^D~ and B^ D^D**~ decays could be 



measured in a time- dependent Dalitz plot analysis of B^ D^D 7r° 20|, |4l|, or the inter- 



ference between CP-even and CP-odd helicity states could be probed in a time- dependent 



analysis of P° D*+D*~ [271, l29|. Both these techniques require the reconstruction of a 



high multiplicity final state as well as precise understanding of potential misreconstruction e 



fects. Another interesting possibility is provided by the vector- vector decay P° J/ipp^ 29| . 



though due to the large natural width of the p, a time- dependent analysis of B^ J/ iI!'k^ti 
may be necessary to incorporate correctly all interference effects 42|]. This would be a highly 
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challenging analysis, though it could potentially be studied at LHCb. 

Another interesting possibility that has been^ proposed involves a study of helicity am- 
plitudes in the dibaryon decay 



However, this decay has not yet been 



observed [431], meaning that it will be difficult to accrue sufficient statistics for a precise 
analysis. 
b qqs 

Since measurements of mixing- induced CP violation phenomena in 6 — > qqs transitions 
provide one of the most interesting approaches to search for effects of physics beyond the 



Standard Model 



|i6l,Q,Q,[~^ 



,|19|], it is clearly important to be able to probe cos(2/3) in these 
transitions. Recently, this has been achieved using time-dependent Dalitz plot analyses of 
KgK~^K~ (containing contributions from (pKg and foKg among others) 44| and 



45, 



The 



Kgir^n^ (containing contributions from p^Kg and foKg among others) 
approach that has been adopted is to obtain values of /3, rather than cos (2/3) and sin (2/3) 
separately. The current results indicate that values of P closer to the Standard Model 
solution than those with cos(2/?) < are preferred, but much more precise results are needed. 
Although LHCb is expected to make some improvement on the current measurements, these 
channels provide one of the motivations for a very high-luminosity electron-positron flavour 



factory [12 



Q,Q,y. 



In Table H] we summarise the current status of experimental measurements of cos(2/5). We 



have not included results from time-dependent Dalitz plot analyses of B^ KgK~^K 44] 



and B 







45 



[46 1 where the results have been presented in a different format. 
Discounting the results in — > D*^D*"Kg, which suffer from a large theoretical uncer- 
tainty, we see that no measurement has a precision better than about 0.50, and therefore 
additional approaches are very welcome. 

In Table HTl we summarise the existing measurements of sin(2/5) in b cud transitions. 
Additional approaches that can improve the precision beyond that achieved in B'^ D^ph^ 
would be very helpful to test the Standard Model prediction that the value of sin (2/3) 
measured should be the same as that obtained from B^ J / tpK^ decays. 
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TABLE I: Summary of measurements of cos(2/3). For more details and world averages, see 
Ref. ^. For all quoted results, the first uncertainty is statistical and the second systematic. 
Note that the parameter measured in analyses of — > D*~^ D*~ K'^ marked by (*) is not cos(2/3) 
but (2Js2/Jo) cos(2/3), where the combination of hadronic parameters {2Js2/ Jo) is expected to be 
positive. In the BaBar results on DhP with D — > ii'^vr+vr", the third uncertainty is due to 

the D decay model; Belle include these effects together with other systematic uncertainties. The 
symbol denotes a light neutral meson such as a vr*^. 



Experiment cos(2/5) 
S° ^ J /ipK* (892) with K*(892) ^ K'^tt^ 
BaBar [26] 3.32 toll ± 0-27 

Belle [32] 0.56 ± 0.79 ± 0.11 

BaBar [35] (*) 0.38 ± 0.24 ± 0.05 
Belle [36] (*) -0.23 Ig'j^ ± 0.13 
5° ^ Dh^ with D K^TT+Ti^ 
BaBar [40] 0.42 ± 0.49 ± 0.09 ± 0.13 
Belle [3^ 1^-^^^^ 



TABLE IL Summary of measurements of sin(2/3) in 6 — > cud transitions. For more details and 

n 

world averages, see Ref. [5[. For all quoted results, the first uncertainty is statistical and the second 
systematic. In the BaBar results on B^ Dh? with D —y ^C^vr+vr", the third uncertainty is due 
to the D decay model; Belle include these effects together with other systematic uncertainties. 

Experiment sin(2/3) 

B^ ^ Z)g>o 
BaBar [4r] 0.56 ± 0.23 ± 0.05 
S° ^ Dh^ with D K'^gTT+TT- 
BaBar [40] 0.29 ± 0.34 ± 0.03 ± 0.05 
Belle [39] 0.78 ± 0.44 ± 0.22 
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III. METHOD 



Detailed descriptions of the method to extract CP violating phases from B — > Dh^ 



(where h'^ is a light neutral meson such as vr^, r/, etc.) can be found elsewhere 38|, |48|, |49 |. 



(Similar discussions where is a K% meson, relevant for the extraction of the angle 7, can 



also be found in the literature 



50 



5l|.) Here we provide only an outline of the method. 



To simplify the discussion, we initially treat B decays to charm as being flavour-specific - 
ie. we neglect h — * ucd amplitudes, that are suppressed by a factor of approximately 0.02 



compared to the favoured h 
amplitudes at a later stage. 



cud amplitudes 



52| . We consider the effects of suppressed 



Consider the amplitude for a B^ decay to a point in the D^-k^ti Dalitz plot, described 



by the coordinates mi = rn?{DTi^) and rn?_ = im?{DTi ). We define the amplitude as 



(1) 



where we have used the isobar formalism to express the amplitude as a sum of contributions 
from interfering resonances. The complex coefficients Cj describe the magnitude and phase of 
the contribution from each resonance i, while the strong dynamics (lineshapes and angular 
distributions) are contained within the Fi functions. The sum over i will include excited 
D mesons [eg. -D2~) ^"^^ tt^tt^ resonances {eg. p°). Feynman diagrams representing these 
amplitudes are shown in Fig. [H 





FIG. 1: Feynman diagrams for (left) a colour-favoured h cud transition resulting in a — >■ 
D2~tt'^ decay, and (right) a colour-suppressed b cud transition resulting in a i?*^ — > decay. 

The amplitude for the decay to D^7i^7i~ can be written similarly, 

A(5° ^ D°7r+7r~) = A{ml, ml) = ^ Qi^i(m^, ) , (2) 
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where the index i is understood to run over the same set of resonances, except that the 
charge of the excited D mesons will be opposite compared to the decay case. Since F 
contains strong dynamics only, we can assert Fi{m'^,m'^) = Fi{m^_,rn?^) and, neglecting 
direct CP in B decay, we have Ci = Ci both for the excited D mesons and for the vr+vr" 
resonances (which can only be CP-even). 

Following the usual formalism (see, for example, 

a 0.3) 

, we write the time-dependent 

decay rate of a -B meson that is known to have specific flavour content (being either B^ or 
B^ at time At = 0) to a particular final state / or its CP conjugate / as 

r(S°hys ^ f{^t)) e-l^*l/^so (1 _ Sf sin(AmAt) + Cf cos(AmAt)) , (3) 
r(^phys ^ /(^^)) e-l^*l/"BO (1 + Sf sin(AmAt) - Cf cos(AmAt)) , (4) 

where Sf = 2Im(A/)/(l + |A^|), C/ = (1 - |A}|)/(1 + |A}|) and A/ = In the above, Tbo 
is the average lifetime of the neutral B meson. Am is the mass difference between the two 
eigenstates of the B^-B^ system which are described in terms of the flavour specific states 
by the mixing parameters q and p as Bk^h)^ = p \B^) + (— ) q B^'j. (We have neglected the 
lifetime difference and assumed CPT invar iance.) 

In the case at hand, / represents a point in the Dti^tt" Dalitz plot. If we consider D 
decays to CP eigenstates (and neglect direct CP violation in the D system), then we have 

^^'"+''"-^"/^A(m^,m2_)' 

where ?7d is the CP eigenvalue of the Dcp state. With the further substitution \q/p\ = 1, 
arg(g/p) = —2/3, both good approximations in the Standard Model, we can write 

2 2n 2im{e-^^^vnA*A) _ 2 Im (e-^^^ E^ c*F,{ml, nilY Ej CjFjjml, ml)) 
b{m^,m_)- |^|.+|^|^ - \j:.c,F,{ml,ml)\^ + \E^c,F,iml,ml)\^ ' 





A 


2 


EiQF,(m^,m2)|2- 


j:icMml,ml)\^ 


\A\^ + 


A 


T — " 


EiCiF,{ml,ml)\^ + 


E,cA{ml,ml)\^ 



C(m+,m„)= —rjj^ = F^^2' 2M2 , ^^^2' 2M2 " (7) 



The numerator of the expression for S{m'^,m1) can be written as 
2?7£) (|cos(2/5)Im(A*A) — sin(2/3)Re(A*A)^ , making explicit the dependence of mixing- 
induced CP violation on both sin(2/5) and cos(2/3). We note that one can choose to fit for 
sin(2/3) and cos(2/3) independently, or alternatively one can fit directly for /3. Although 
the latter appears attractive, the expressions above make clear that in the former case 
both parameters appear as coefficients of physically observable functions of the amplitudes, 
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and thus one might expect somewhat better statistical behaviour for these observables, in 
particular in regions close to physical boundaries. 

The sensitivity to cos(2/5) is proportional to Im(y4*y4), and therefore depends strongly 
on interference in the Dalitz plot between resonances with non-trivial phase differences. 
To illustrate this point, we show in Fig. [2] the values of Im(74*74) and Ke{A*A) across the 
DcpTc^TT^ Dalitz plot, calculated using our nominal model as described in the next 
section. The Dalitz plot is drawn as = "^dtt- ^+ = ''^'dtt+- emphasise that the 
Dalitz plot model can be experimentally determined from the higher statistics D^7t^7t~ 
sample with flavour-specific -D° decays. Consequently, the model uncertainty on the obtained 
value of cos(2/5) should be controllable. 




5 10 15 20 25 5 10 15 20 25 



ml (GeV V) ml (GeV^/c") 

FIG. 2: Dalitz plot distributions of (left) |lm(A*^)| and (right) |Re(A*74)|, which govern the 
sensitivity to cos(2/3) and sin(2/?), respectively, using our nominal — > Dcpn'^TT" decay model. 
Note that the z-axis is shown on a log scale. 

We conclude this section with a brief discussion of some of the approximations that have 
been made in the above formalism. We have assumed that the B^-B^ system contains no 
CP violation in mixing {\q/p\ = 1), negligible lifetime differences (AF^ = 0) and no CPT 
violation. Furthermore, we have assumed no direct CP violation in D meson decays or in B 
decays to Dtt^ti^. All of these are valid approximations in the Standard Model, and have 
been experimentally tested to good precision. 

We have also until now neglected contributions from the suppressed b ucd amplitudes, 
illustrated in Fig. [31 The relative weak phase between h cud and h ucd amplitudes is 
given by the CKM Unitarity Triangle angle 7. If the contribution of the suppressed decays 
is significant, it leads to some interesting phenomenology, including the potential to measure 
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sin(2/3 + 7) from resonant amplitudes such as D'^'^tt^ , and to measure 7 from rates and direct 
CP violation effects in modes such as Both measurements would, however, benefit from 
the larger statistics that are available by reconstructing the D meson in a flavour-specific 
decay mode and are anyway unlikely to be competitive with similar measurements using 
/}{*)±7j-T Qjid DK final states, respectively. This serves to illustrate how flavour-specific 
D decay modes can be used to control model uncertainties, as well as other experimental 
systematic uncertainties, in the analysis. A further corollary is that, unless the suppressed 
amplitudes are accounted for, there can be small biases on the extracted values of 2(3 that are 
measured from Dcptt^tt" , with the biases opposite in sign for D mesons reconstructed 



in CP-even and CP-odd final states 49|]. The suppressed amplitudes therefore introduce 
some model dependence into the results, which can be tested by adding suppressed ( "wrong- 
sign") D2 resonances (for example) into the model, by relaxing the constraint Cj = Cj, and 
by checking the consistency of results of independent fits to the samples with CP-even and 
CP-odd D decays. 





d 

P" 



FIG. 3: Feynman diagrams for Cabibbo-suppressed transitions: (left) b — > ucd resulting in 
L'2^7r~, and (right) b — > ■ucd resulting in — > D^p^. 
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IV. FEASIBILITY STUDY 



A. Nominal Model 



We form a nominal D^ti^-k Dalitz plot model based on results from the Belle 

Collaboration 22|. The most prominent contributions to the Dalitz plot are found to be 
5° ^ 01-11+ , Dif D^TX-- 5° ^ , Dl~ ^ D^TX-- fiO ^ Dl~ii+ , D*~ D^tt^; 

D^p^ , 71+ and D^f2 , tt+tt^. The parameters of these resonances 

are taken from the Particle Data Group (PDG) 53|] and are summarised in Table lllll The 
particle denoted Dl~ is a virtual D*~~ meson. The D*" itself is too long-lived to cause 
interference and, following Belle 22|], we veto the region dominated by this contribution by 
excluding the invariant mass range 2.00GeV/c^ < mj:,^^ < 2.02GeV/c^ from the analysis. 



TABLE III: Parameters of the resonances used in our nominal model. These vahies are taken from 
the PDG [531. 



Resonance 


Mass (MeV/c^) 


Width (MeV) 


n*~ 


2461.1 ± 1.6 


43 ±4 




2352 ± 50 


261 ± 50 


D*~ 


2010.27 ±0.17 


0.096 ± 0.022 




775.49 ± 0.34 


149.4 ± 1.0 


f2 


1275.1 ± 1.2 


185.0 l|j 



The model described above is sufficient for the purposes of our feasibility study, although 
the D^-K+Tx- decay amplitude may include additional resonant or nonresonant terms. Belle 
found a possible contribution from B^ Z)''/o(600) (the /o(600) is sometimes known as 
the a meson). Ultimately the decay model and its uncertainty should be determined from 
n+n with flavour-specific decays as a part of the analysis. 

B. Event Generation 

We estimate the number of B^ Dcptx+ti- events that we expect in the final BaBar 
dataset as follows. Belle found 2909 ± 115 5*^ — D^tt+tt- , K+tt~ events in a data 
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sample of 388 x 10^ BB pairs. We assume similar selection efficiency to that in the Belle 
analysis, and scale according to the size of the final BaBar dataset: 467 x 10® BB pairs. 
We calculate the expected numbers of events in CP-even {K^K~ and ti^-k') and CP-odd 
(i^gvr^ and K^g^) events using scaling factors calculated from event yields in a recent analysis 
of P+ DK+ from BaBar [sj] . These factors are (8.6 ±0.2)% for K+R- and (3.1 ±0.1)% 
for TT^ir" , giving a total (11.7 ± 0.2)% for CP-even; correspondingly (8.9 ± 0.2)% for Kgn'^ 
and (3.3 ± 0.1)% for K'^oo sum to a total (12.2 ± 0.2)% for CP-odd. Taking all these factors 
into account, we expect approximately 410 CP-even events and 430 CP-odd events in the 
final BaBar data sample. 

We use the Laura++ package to generate and fit events. This package has been largely 
developed by the authors and used in several recent analyses published by the BaBar Col- 



laboration 
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56 



57, 



58 



59|. We generate events using q parameters based on results of 



the Belle analysis of P° D^ti^tt^ |22| and = Cj, using a convention in which all Pj 
functions are normalised to unity when integrated over the Dalitz plot. We assume a recon- 
struction efficiency that does not vary across the Dalitz plot and neglect misreconstruction 
effects. Initially, we also neglect backgrounds from other B decays or other processes, but 
we include these events in the study at a later stage, discussed below. We simulate effects 
due to At resolution and misidentification of the flavour of the tagging B meson using stan- 
dard resolution functions and parameters from BaBar j6o|. In the fit we float the real and 
imaginary parts of all q except those for the which are fixed as reference parameters. 
We also float both cos(0mix) and sin(0mix)j where (pmix = 2/3 in the Standard Model. The 
generated distribution of events in the Dalitz plot for ten times the expected statistics can 
be seen in Fig. HI Projections onto m^, and 7717^+^- can be seen in Fig. [51 



C. Results 



The distributions of the fitted results for cos((/)mix) and sin(</)mix) from 500 pseudo- 
experiments generated as described in the previous subsection are shown in Fig. [61 We 
are clearly able to determine both cos(0mix) and sin^cp^i^), with spreads of the distributions 
of about 0.50 and 0.18 respectively. (For comparison, the means of the distributions of the 
uncertainties on these parameters that result from the fits are found to be 0.43 and 0.17 
respectively, conffiming the approximately Gaussian nature of these parameters that is also 
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ml (GeV^/c^) 



FIG. 4: Dalitz plot distribution of the generated events using our nominal model. The statistics 
shown here are ten times that expected in the full BaBar dataset. 
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0.2 0.4 0.6 0. 



'"D^.aGeV/c^) 



m_ (GeV/c^) 



FIG. 5: Projections onto (left) nij^^p^± and (right) m^+^- of the generated events using our 
nominal model. The statistics shown here are ten times that expected in the full BaBar dataset. 
To avoid artefacts due to reflections, in the mj^^p^i projection we require m^+^- > 2.0 GeV/c^; 
in the m^+^- projection we require rnjy^p^± > 2.75 GeV/c^ (both combinations). Structures due 
to (left) ^-iid (right) and /2 resonances are clearly apparent. 

evident in Fig. [61) Small biases in the fit results are found to disappear when the number 
of events per experiment is increased. All fitted q parameters are found to be similarly 
consistent with the input values. 

Further to illustrate the method, we show in Fig. [7] the At asymmetry for events in 
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FIG. 6: Distributions of the fitted values of (left) cos((/)niix) and (right) 

sin((^i^ix) from the pure 

signal study. The resolutions are 0.50 it 0.02 and 0.18 it 0.01, respectively. The red arrows indicate 
the generated values, which are based on the value 2(3 = 21.7°. 

the region of the p° resonance in the distribution. The asymmetry is between events 

where the other B meson has been identified ( "tagged" ) as a 5° or a by the charge of the 
lepton produced through its semileptonic decay. If this region contained B^ DqpP^ decays 
alone, the asymmetry would give sin(2/5) sin(AmAt), smeared by experimental effects (At 
resolution and misidentification of the flavour of the tagging B meson). The plot contains 
only events containing a D meson reconstructed in a CP-even decay mode and has 500 
times the statistics expected from the final BaBar dataset. 




-10 -8 -6 -4 -2 2 4 6 8 10 



A t (ps) 

FIG. 7: Asymmetry between and tags as a function of At in the region of the resonance 
(0.70GeV/c^ < 771^+^- < 0.85GeV/c^). Only events in the best tagging category are shown. The 
statistics shown here are 500 times that expected in the full BaBar dataset. 

To ensure that our results do not depend on the true value of 0niix we repeat the process 
with a number of different input values. As shown in Fig. [HI we determine the correct values 
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of cos(0mix) and sm(0rnix) in all cases. We also find that uncertainties on these parameters 
do not depend strongly on the input values - a variation of 4% (18%) is found in the mean 
uncertainty on the fitted value of cos(0mix) {sm{(f)jnix)) , with the uncertainty being largest 
when 0mix is smallest. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Generated cos(0 ) Generated sin(0 ) 

mix mix 



FIG. 8: Fitted values as a function of the generated values of (left) cos((/)niix) and (right) sin((/)mix)- 
In the latter plot the error bars are too small to be seen. 



D. Experimental Complications 

In a true experimental environment there will be events from background processes that 
will have a specific structure in the Dalitz plot and also in At. The presence of these 
events will complicate the analysis and could degrade the sensitivity to the parameters of 



interest. In order to attempt to estimate the sea 



including background events. The Belle analysis 22] indicates that the level of background 



e of these effects we repeat the study 



is approximately the same as that of signal once they apply a selection on the discriminating 



kinematic variables tties and AE (for definitions of these variables, see for example |58| ) . 
As such we include the same number of background events as signal in our samples. We 
take the distribution of the background events in the Dalitz plot from the Belle paper and 
use a delta function for the true At distribution. In order to provide further discrimination 
between signal and background we also include rriEs and AE in the fit. We use Gaussian 
shapes to describe the signal distribution of both of these variables, whilst for background 



we use the ARGUS shape 



611 ] and a linear function, respectively. 



The results of this extended study are shown in Fig. O The fitted values of cos(0mix) and 
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sin(0jnix) are still unbiased and the uncertainties are largely unchanged from the pure signal 



case. 




FIG. 9: Distributions of the fitted values of (left) cos(0mix) and (right) sin((/>niix) from the signal 
and background study. The resolutions are 0.54ib0.02 and 0.16ib0.01, respectively. The red arrows 
indicate the generated values, which are based on the value 2/? = 21.7°. 

There are a number of additional potential experimental complications that we have 
not simulated. Effects due to smearing of the reconstructed Dalitz plot position should be 
negligible since our Dalitz plot model does not contain any very narrow resonances. More 
significant misreconstruction effects can occur when one of the particles from the signal side 
is exchanged with a particle from the decay of the other B in the event. This effect tends 
to occur most frequently near the corners of the Dalitz plot. Since we veto the regions in 
the corners of the Dalitz plot that are dominated by decays of the D* meson, we do not 
expect this to cause any serious difficulty. Finally, we have assumed that the reconstruction 
efficiency does not vary across the Dtt^tt^ phase space, but in a real experiment it is likely 
that the corners of the Dalitz plot have a lower efficiency. This may lead to a small reduction 
in the precision of the cos(2/5) measurement, since these are the regions that are most 
sensitive to this parameter, as shown in Fig. [21 However, in an experimental analysis the 
efficiency would be measured using detailed Monte Carlo simulation together with data 
control samples, and therefore we do not expect any bias on the results. 



E. Prospects at LHCb 

The potential to utilize this method at LHCb has not been explicitly studied, but some 
useful extrapolations can be made. We consider only the CP-even D meson decay to K^K~, 
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which is well-suited for study at LHCb since it has a final state consisting of only charged 
tracks including two kaons. (The use of particle identification information from LHCb's 
ring imaging Cherenkov detectors being essential to reduce combinatoric background in the 
hadronic environment 6^.) 

An estimation of the likely yield can be made by comparison with the decay D^n^, 
— s> K^7r~iT^ . In this channel, which has a product branching fraction of (2.47 ± 0.13) x 
10- LHCb expects to trtoer and reeonstruct 1,34 x 10» events in 2 fb- of data (o„e 
nominal year of data taking) 63|]. For Dir^n^ , D K~^K^ the product branching 

fraction is (3.3 ±0.4) x 10^^ so that if the trigger and reconstruction efficiencies are the same, 
one ex pec ts approximately 18,000 events. Taking an effective tagging efficiency at LHCb 
of 5% 64], the equivalent number of perfectly tagged events in one year of LHCb data is 
approximately 900, which compares very well to the yields from the final BaBar dataset 
(where the effective tagging efficiency is about 30% |60|]). Although it may be necessary 
to apply less efficient selection criteria to suppress background, this provides an indicative 
measure of the potential at LHCb. With 2 fb^^ of data, it will be possible to achieve 
a precision better than any previous measurements. With the complete LHCb data set, it 
should be possible to measure cos(2/3) and sin(2/5) with precisions of 0.12-0.17 and 0.03-0.05 
respectively. 

The prospects for LHCb to make precise measurements of cos(2/9) and sin(2/3) using 
B'^ Dn^Ti" , D therefore look rather good. However, any firm conclusion on 

this point requires a detailed study including proper simulation of detector effects as well as 
consideration of the Dalitz plot model including effects of the suppressed amplitudes. We 
leave such studies to further work. 



18 



V. CONCLUSION 

We have presented the results of a feasibihty study of a method to measure cos(2/3) using 
a time-dependent Dahtz plot analysis of Dtt^tt^ decays, where the neutral D meson is 

reconstructed in decays to CP eigenstates. We estimate that, with the final BaBar dataset, 
cos(2/5) can be measured with a precision of ~ 0.50, making this approach competitive with, 
or superior to, all other methods that have been attempted to date. Furthermore, sin(2/3) 
can be measured to within ~ 0.16, which is more precise than any existing measurement 
using b — > cud transitions. 

We have argued that uncertainties relating to the composition of the Dtt+tt" 
Dalitz plot can be tamed using the much larger data samples that are available when the D 
is reconstructed in a flavour-specific hadronic decay mode such as K^n^. Therefore, 

the sensitivity is limited by statistics only, and more precise measurements will be possible 
at electron-positron colliders with higher luminosity. Moreoever, since final states containing 
only charged particles can be used, there is great potential for this analysis to be employed 
at LHCb, where use of the CP-even decay D K^K~ looks particularly promising. 

Application of this method at future experiments will allow a determination of cos(2/?) 
that would definitively establish its sign, hence resolving the ambiguity on 2(3. Moreover, 
measurements of sin(2/5) in 6 ^ cud transitions can be made with precision comparable 
to that obtained by the B factories in 6 ^ ccs transitions (5° J/ipK^). Discrepancies 
between these values would be unambiguous signs of new physics. 

This work is inspired by the luminosity recorded by the B factories, and we are par- 
ticularly grateful to our colleagues from PEP-II and BaBar. We would like to thank Alex 
Bondar, Paul Harrison, Gagan Mohanty, Guy Wilkinson and Jure Zupan for reading the 
manuscript and making useful suggestions. This work is supported by the Science and 
Technology Facilities Council (United Kingdom). 
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